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The Auxiliary Field Diffusion Monte Carlo method has been applied to simulate droplets of 7 
and 8 neutrons. Results for realistic nucleon-nucleon interactions, which include tensor, spin-orbit 
and three-body forces, plus a standard one-body confining potential, have been compared with 
analogous calculations obtained with Green's Function Monte Carlo methods. We have studied 
the dependence of the binding energy, the one-body density and the spin-orbit splittings of on 
the depth of the confining potential. The results obtained show an overall agreement between the 
two quantum Monte Carlo methods, although there persist differences in the evaluation of spin- 
orbit forces, as previously indicated by bulk neutron matter calculations. Energy density functional 
models, largely used in astrophysical applications, seem to provide results significantly different from 
those of quantum simulations. Given its scaling behavior in the number of nucleons, the Auxiliary 
Field Diffusion Monte Carlo method seems to be one of the best candidate to perform ab initio 
calculations on neutron rich nuclei. 



I. INTRODUCTION 

Neutron rich matter is a subject of fundamental inter- 
est in nuclear astrophysics. Important phenomena, hke 
the structure and evolution of compact stellar objects, 
the r-process in nucleosynthesis or the mechanism of su- 
pernovae explosion cannot be understood without a deep 
knowledge of the properties of such matter^, 13,11,0)0 ■ 

In recent years the relevance of nucleon-nucleon (NN) 
correlations on important quantities, like the equation of 
state, density and spin responses or neutrino opacity, has 
been firmly recognized^ 0,11, 0^3- Mean field models 
where neutron matter is treated as a weakly interacting 
system are insufficient to describe these correlations. 

Dense and cold nuclear matter cannot be reproduced 
in our terrestrial laboratories. Therefore, it must be stud- 
ied theoretically, and one needs to make use of the full 
machinery of modern many-body theories. 

Exotic nuclei with large neutron excess, far from the 
stability valley, including isotopes of nitrogen, oxygen 
and fluorine, have been produced in heavy ion reactions 
with radioactive beams, and subsequently studied spec- 
troscopically. Experimental data on their binding ener- 
gies and neutron removal energies provide a very impor- 
tant body of information on the structure of neutron rich 
matter and on the NN interaction, at least in the low den- 
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sity regime. It is well known that tensor, spin-orbit and 
three-body forces play a fundamental role. Small modi- 
fications of these forces may lead to large effects on the 
properties of dense and cold hadronic matter. 

Neutron droplets were originally introduced as a home- 
work problem to compare ab initio calculations with en- 
ergy density functionals commonly used in astrophysical 
investigations. They also provide clean benchmarks for 
modern many-body theories, being simpler systems than 
either nuclei or nuclear matter. 

Closed shell drops made of eight neutrons, and open 
shell ones with seven or six neutrons have been studied 
in the last decade, by using quantum Monte Carlo meth- 
ods, like Variational Monte Carlo and Green's function 
Monte Carlo (GFMC)^3l- These quantum simulations 
carried out with a realistic bare NN interactions deviate 
significantly from energy density functional resultsfl^. 

Neutron droplet models can also serve to study neu- 
tron rich nuclei, devising one-body effective potentials 
to describe the interaction of the halo neutrons with an 
inert core. For instance, the oxygen isotopes can be mod- 
eled by neutron droplets with up to 10 neutrons out of 
the ^^O core, and fluorine isotopes with droplets made of 
one proton and up to 12 neutrons in the halo 13J. 

The recent developments made in Quantum Monte 
Carlo methods, allow for ab initio calculations of 
medium-heavy nucleonic systems, at an unprecedented 
accuracy. GFMC calculations of the binding energy 
of nucleonic systems with A = lA have been recently 
been performed 14j. Auxiliary Field Diffusion Monte 
Carlo (AFDMC) requires order operations per time 
step and this polynomial scaling allows the study of 
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many more nucleons. For example, neutron matter 
calculations [isj have used 114 neutrons in a periodic box 
to calculate the equation of state. In both methods the 
spin-isospin dependence of the nucleonic interaction has 
been handled without any approximation. 

In particular, the AFDMC method looks very promis- 
ing to perform quantum simulations of large nuclei and 
nuclear matter 15, 16] . The main reason is because it 
samples, rather than performing a complete sum of the 
spin states. The sampling is not done directly on the spin 
states, but auxiliary fields of the Hubbard-Stratonovich 
type are introduced to linearize the spin dependence of 
the NN potential operator. The method consists of a 
Monte Carlo sampling of the auxiliary fields, and then 
propagating the spin variables at the sampled values by 
means of a rotation of each particle's spin spinor. 

In this paper we apply the AFDMC method to study 
neutron drops of 7,8 neutrons, kept bound by an effec- 
tive one-body potential of the type given in ref. 12]. The 
main goal is to test the newly developed AFDMC ver- 
sus other quantum Monte Carlo methods, like GFMC 
or Cluster variational Monte Carlo. Reasonably good 
agreement is obtained for the binding energy of *n(0''') 
and ^n(3/2~), whereas a larger discrepancy exists on the 
spin-orbit splitting of ^n. 

Our calculations represent a preliminary step towards 
the study of neutron rich nuclei. 

The plan of the paper is the following. In the next 
section we describe the hamiltonian used in our calcula- 
tions. A brief outline of the AFDMC method is given in 
Section III. A reformulation of the three-body potentials 
of the Illinois type is given which allows for efficient nu- 
merical calculation. Results are presented and discussed 
in Section IV. The last Section is devoted to conclusions 
and future perspectives. 



II. HAMILTONIAN 

The ground state properties of the neutron droplets are 
computed starting from a nonrelativistic Hamiltonian of 
the following form: 



H 



T + Vi + V2 + V:i 



^ 2r 



11% 



(1) 

i<j<k 



The one-body potential Vext stabilizes the drop, other- 
wise unbound. We take it of the same form as in ref. 
[l3l. It consists of a Wood-Saxon well of the form 



Vextir) 



Vo 



1 + e^"^-^^/"- 



(2) 



where R = 3fm, and a = 0.65fm. The parameter Vo 
determines the depth of the well. We have made calcu- 
lations for four different values of Vq = 20, 25, 30 and 35 



MeV, where the first value has been used in the GFMC 
and variational Monte Carlo calculations of ref.flTl IT^ 
and more recently in ref.^3]- The larger values of Vq 
have been considered to test the reliability of AFDMC 
for more bound systems, with a large peak density. To 
increase the density further R must be decreased. 

The two-body NN interaction considered belongs to 
the Urbana-Argonne vi potentials : 



vi 



(3) 



truncated to include only the following 8 operators 
K)[ll|ll|li: 



0^=^^\i,j) = {l,<j,-<jj,S,,,L,, 



S-.j)® (l,fi-fj) (4) 



3(7; 



(T,- • (To is the 



where the operator 5y 
tensor operator and Lij ~ —ifirij x (V; — Vj)/2 and 
Sij — h(ai + (Jj)/2 are the relative angular momentum 
and the total spin for the pair ij. For neutrons fi-Tj — 1, 
and we are left with an isoscalar potential. 

The Vg potential is a simplified version of the uig po- 
tential, having the same isoscalar parts of uig in all S 
and P waves, as well as in the ^Di channel and its cou- 
pling to the ^Si. It is only semirealistic because it does 
not fit the Nijmegen N-N data ^3] at a confidence level 
of x^/^data ~ 1, as Wis does. However, the difference 
between wig and v'^ is rather small for densities smaller 
or of the order of the nuclear matter equilibrium density 
Pq ~ O.lGfm"'^, and it can be safely added perturbatively. 

The v'g potential should be considered as a realistic 
homework potential, and it has been used in a number 
of calculations on light nuclei llyl Il9l. symmetric nuclear 
matter 21], neutron matter [l5Ll2lIl22 ] and spin polarized 
neutron matter^;|]. 

In order to estimate the contribution of the spin-orbit 
force, we have also made calculations with a two-body 
potential obtained from dropping the spin-orbit terms. 
We denote this potential as Wg, but it should not be con- 
fused with the v'q potential given in ref. [2^ which gives 
a correct binding energy of the deuteron. 

The three-body interaction considered encompasses 
the form of both the Urbana and Illinois 3-body 
potentials^j. Results will be given for the Urbana IX 
potential. The Wg or the v'q two-body potentials plus the 
Urbana IX interaction are denoted as AU8' and AU6' 
respectively. 



III. AUXILIARY FIELD DIFFUSION MONTE 
CARLO METHOD 

The Auxiliary Field Diffusion Monte Carlo is an 
extension to DMC to deal with spin dependent hamilto- 
nians. The quadratic dependence of these hamiltonians 
on the spin operators is taken care of by sampling auxil- 
iary variables, which serve to linearize such dependence 
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through Hubbard-Stratonovich transformations. A de- 
tailed discussion of the method can be found in Ref. 
[T5I I25I I . Here we hmit ourselves to briefly outlining the 
method. 

The Vq two-body potential can be separated into a 
spin-independent and a spin-dependent part: 



with 



V = 



ySI _^ ySD ^ 

0'iaAia:jf30'jfj , 



E 



where the elements of the matrix A are given by the 
proper combinations of the components Vp in Eq. ||2J). 
Latin indices, like i and j, are used to for particles, while 
the greek ones, like a and /3, refer to the Cartesian com- 
ponents of the operators. We use the summation conven- 
tion that all repeated greek indices are summed from 1 
to 3. 

Because Aia-.ip = the 3A^ by 3A^ matrix A has real 
eigenvalues and eigenvectors, defined by: 



(5) 



The spin-dependent potential can therefore be written 
in terms of such eigenvalues and eigenvectors in the fol- 
lowing form: 



J/3 



If one defines new N-body spin operators 

On = ^CTiaV'^ , 



the spin-dependent potential becomes 

2 



V 



, 37V 



(6) 



(7) 



(8) 



In the short-time limit we can decompose the imaginary 
time propagator of the diffusion process, which projects 
the ground state out of a trial wavefunction in the fol- 
lowing way: 



-HAt 



-TAt ^-VcAt -V^" At 



(9) 



where Vc = ^ Vextifi) + V^^ is the spin independent 
part of the interaction. The propagation accounting for 
the kinetic and Vc operators gives rise to the usual drift- 
diffusion scheme of DMC. The spin-dependent two-body 
potential part e^^ '^'^ is handled by making use of the 
following Hubbard-Stratonovich transformation 

1 



e ^ 



-foo 



dxn e 2 



, (10) 



n 



--A„O^Ar 
e ^ 



(11) 



where the commutators amongst the 0„ are neglected, 
which requires to keep the time step Ar small enough. 

In Eq. (|1U|) the quadratic dependence on the spin oper- 
ators is transformed into a linear expression which cor- 
responds to a rotation in the spin space. For each eigen- 
value A„ a value of Xn is sampled, and the current spinor 
value for each particle is multiplied by the set of matrices 
given by the transformation in Ea. (|10|) . 

The spin-orbit and three-body potentials can be 
treated within the same scheme. It is important to no- 
tice that while the spin-orbit potential is already linear 
in the spin operator, it is necessary to eliminate spurious 
terms from the simple linearization of the propagator in 
order to take into account corrections at order Ar. This 
leads to additional two- plus three-body counter terms, 
which can be treated as additional interaction terms 1151. 



A. Three— body potential 

We give in this section a reformulation of the Urbana 
and Illinois 3-body potentials 17] for neutrons, which pro- 
vides an efficient and simple way to program them. They 
can be written so that the 0{N^) parts are done with ma- 
trix multiplies. All of these potentials can be written in 
the form 



Vijk — ^277 '-^ijk + ^2t: ^ijk 

The 1-pion exchange amplitude in Argonne wis is 
1 f 



tNN 



where 



X 



op 



3 4^ '^^N. ' 



(12) 



(13) 



+ Y{m^rij)ai ■ gj , (14) 
and the functions Y{x) and T(x) are given by 



Y{x) 



nx) - {^^+^- + l]Y{x)^T{r) 
^Y{r) = ^T{r) = 1 - e 



(15) 



The operator ^ - j X°j' has the same algebraic struc- 
ture as the spin dependent part of the two-body potential 
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V^°. Therefore, it can be expressed in a similar way, in 
terms of a 'UN by 3A'^ matrix Xi^-jp as 



(16) 



with Xia-^ifi — 0. Notice that this matrix is symmet- 
ric under cartesian component interchange a <-*■ /? and 
under particle label interchange i <-» j and is also fully 
symmetric X^p-ia = X^a-jp- 

It is convenient to calculate the square of the X matrix 



(17) 



The Fujita-Miyazawa form for O^^il^^ and neutrons 



becomes 



E 0--^'^=4^a.„a,,X,i,,. (18) 

i<j<k i<j 

The spin independent O^^. terms can be written using 
just pair sums as 



(19) 



i<j<k 



The Tucson S wave component can be written in terms 



of 



Zix) = -[Yix)~Tix) 
= 0. 



(20) 



The two-pion, S wave component of the potential be- 
comes 

i<j<k i<j k 

Note the inner sum makes this a matrix multiply. This 
and the X'^ calculation above are the only order op- 
erations needed. 

The 3-pion terms have a central part and a spin de- 
pendent part. 



i<j<k 

The central part is 



_ 400 ^ 2 y 



(22) 



(23) 



which is 200/9 times the trace of X^. The spin dependent 
part is 

200 ^ A x^2 /r\ a\ 



i<3 



and one can just write out the 4 nonzero terms for each 
combination of and aj . 

As might be expected, the expression used above is 
related to S^f. of ref.[l3 if only the k term of the X^ 
sum is kept, 

'^^jki ~ '^ia(^jl3Xi^.^k(t>Xkip;jKXiS-ji^eajSepf^i^ . (25) 

The / part is 

35^ — 3 y^ S'j^fc -l- "iXia-jpXip.j^Xij-ja ■ (26) 

eye 

Inserting 1 or 3 Levi-Civita s ymb ols instead of 2 will 
produce the A operators of ref. [13 , 

— '^^ia'^j/S^kfiXi^-jj^Xj^-j^i/Xip-j^g 
•^^adjk = —icriaXiyjpXjp-kiiXifj^.^kd^ajS- (27) 

The spin dependent part of the three-body interaction 
can be easily included in the matrix Aia-jp by 



^^27r / , '-^a;ifc'-^. 



13; jk 



+ Xf^.j^XiS-juif-a-yS^I^KW ■ (28) 



B. Trial wave function 



The wave function used as a trial and importance func- 
tion for the DMC algorithm has the following form 



^{R,S)^Fj{R) D{R,S) 



(29) 



where R = (ri, . . . , rV) and S = (si, . . . , s^v). The spin 
assignements Si consist in giving the spinor components, 
namely 



U^ 



U.^ ])+d,\ i) , 



where Ui and di are complex numbers. The Jastrow cor- 
relation operator is given by 



and 



l<3 



D{R,S) = detM 



(30) 



(31) 
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TABLE I: Variational parameters Vq"'' , R""'' and a"'"' of the 
Wood-Saxon potential Vws found in correspondence of the 
various depths Vb of the well of Vext- 



Vo(MeV) 



20 25 30 35 



Vo"'"'(MeV) 20.0 25.0 30.0 35.0 
i?'""'(fm) 1.5 1.6 1.4 1.2 
a"°''(fm) 0.60 0.80 0.80 0.80 



is the Slater determinant where the elements of the Slater 
matrix are one-body spin-space orbitals, 



Mij = {ri,Si 



(32) 



The Jastrow function fj has been taken as the scalar 
component of the FHNC/SOC correlation operator Fij 
which minimizes the energy per particle of neutron mat- 
ter at density p = 0.16 fm~^ 0,^3- The single particle 
orbitals are taken as solutions of the Schroedinger equa- 
tion 



2m 



Vws{r) 



(f) = e(j>{r), 



(33) 



where Vws is a Wood-Saxon well, hke Vext of Eq. (0) 
with parameters chosen to optimize the expectation value 
of the energy 



(*l*> 



(34) 



for each value of Vo considered. Et is calculated with 
variational Monte Carlo. The parameters of the resulting 
Wood-Saxon potential are reported in Table |l| 

In the case of ®n(0"'") the orbitals fill the s and p shells. 
In the case of ^n(3/2^) the missing orbital off the p-shell 
is in the state \l,m,s = 1,1,1/2); for ^n(l/2~) we have 
considered the combination (|1, 1, -l/2>-|l,0, l/2))/v^. 
The above choices give rise to different determinants Dq, 
D3/2 and Di/2 respectively. They are evaluated at the 
current values of the positions and spin assignments of 
the nucleons in the walker \R,S). 

The Jastrow operator Fj misses all the spin-dependent 
correlations. Therefore our trial function is rather poor 
and we do not expect particularly good results for our 
mixed estimates, except for the total energy. Better 
trial functions can be considered, although they require 
a more demanding computational effort. For instance, 
the space-spin orbitals can be modified to include spin 
correlations in the Slater determinant, as recently done 
in neutron matter calculation to take into account the 
spin-orbit correlations in the trial function|27l|. Such im- 
proved trial function provides a better nodal surface, and 
gives a better variance. 



C. The algorithm 

The AFDMC algorithm is implemented as usual, with 
a propagation in imaginary time of a population of walk- 
ers \R, S) according to the propagator in eq. © with the 
standard drift-diffusion procedure. In addition one has 
to sample the a;„ auxiliary variables given in Ec ljlUI) to 
rotate the spinors. After that all the weight factors are 
computed, they are combined to evaluate a new value of 

In order to avoid the fermion sign problem due to the 
antisymmetric character of the wave function, a path con- 
straint is introduced. If the real part of S) is nega- 
tive, the walker is included in the evaluation of the mixed 
and growth energies, but then is dropped from the pop- 
ulation. In general, the importance sampling makes the 
number of dropped walkers small. In our calculations 
here the number of node crossings is of order 1 percent. 



IV. RESULTS 

In Table im we show our AFDMC results for the mixed 
energies, obtained for a ^n{0'^) droplet confined by a 
Wood-Saxon well Vext with a depth value Vb = 20MeV 
in correspondence with the AU8' and AU6' nuclear inter- 
actions. They are compared with the most recent GFMC 
estimates J7J j32| 

One can see that there is an overall agremeent between 
AFDMC and GFMC. The discrepancies on the binding 
energies of ^n{0+) and ''71(3/2-) are within 2%. This 
is very small if one consider that the binding energy re- 
sults from a large cancellation between potential and ki- 
netic energy (the mixed value of the potential energy is 
of — 124.6(3)MeV) and that the variational energy with 
our ^ is only -26.27MeV . 

A larger discrepancy is obtained for the ''n(l/2^) state, 
which makes the AFDMC spin-orbit splitting of "^n about 
half that of GFMC. Such discrepancy may be due to 
differences between AFDMC and GFMC in treating the 
spin-orbit interaction, as alrea dy p ointed out by the neu- 
tron matter calculations of ref . . 

We have attempted to go beyond our path constraint, 
by doing transient estimation. This requires the guiding 
function to have no zeroes where the true ground state is 
non zero. The guiding function that we have used is the 
modulus of the antisymmetric function smoothed with a 
term that decays to zero as function of the distance of 
neutrons from the center of the drop: 



(35) 



where: 



N 



N 



(36) 



a k=l 
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0.000 0.005 0.010 0.015 0.020 0.025 0.030 0.035 0.040 
At (MeV"') 



TABLE II: Ground state AFDMC energies of *n(0+), ^n(i ) 

and ^n(|") droplets for Vb = 20MeV and the AU8' and AU6' 
interactions. The cluster variational Monte Carlo (CVMC) 
and GFMC resuhs of ref.[ll for the AU8' and the full AU18 
(Argonne wig plus Urbana IX) are also reported for compari- 
son. The last column reports the spin-orbit splittings (SOS) 
in MeV of ^n, given by the energy difference between the 
''n(| ) and ''^(i ) states. 



^^71(0+) ) ^n(| ) SOS 



GFMC(AU18) 
CVMC(AU18) 
GFMC(AU8') 
AFDMC(AU8') 



-37.8(1) 
-35.5(1) 
-38.3(1) 
-37.8(2) 



-33.2(1) 
-31.2(1) 
-34.0(1) 
-32.5(2) 



-31.7(1) 1.5(2) 

-29.7(1) 1.5(2) 

-32.4(1) 1.6(2) 

-31.8(1) 0.7(2) 



AFDMC(AU6') -36.82(5) -31.6(2) -30.8(4) 0.8(5) 



FIG. 1: Results of a transient estimation for the ground state 
energy of a **n started from an AFDMC calculation. 



The parameter a has been chosen equal to 0.1 . 
Branching is done on the magnitude of the total weight 
Wk of the k-th walker. The mixed estimate is given by 



O 



mixed 



2^k '^k (■q.^lB.^Sk) 



(37) 



The results obtained for the *n(0"'") drop are shown 
in Fig. 1. As it can be seen, the energy in the imag- 
inary time interval of 0.04 MeV~^ does not show any 
significant decay, and it is compatible within errorbars 
with the estimate of the constrained path energy which 
is —37.8 MeV. The errorbars are in any case rather large, 
which is due to a poor description of spin dependent cor- 
relations in both trial and guiding functions. They can 
be included by using the backflow form of ref. Work 
in this direction is in progress. 

Another interesting comparison is with the results ob- 
tained with various energy density functional models, like 
(i) SkyrmeM 28], (ii) Skyrme 1' 29, 30], (iii) FPSJH and 
(iv) FPS-2lQ, which reproduce the ground state ener- 
gies of stable closed shell nuclei rather accurately. The 
energy density functional models provide a range of val- 
ues for the binding energy of ®ri(0+) which goes from 
-32.1MeVof FPS to -47.4MeV of Skyrme M 12j. More- 
over, the spin-orbit splitting is about 3MeV, too large 
with respect to the quantum Monte Carlo estimates. 

We have repeated the calculations of ground state en- 
ergies and densities in the Skyrme 1' and in the SKM 
models for the ®n(0"'") drop. The interaction parameters 
have been taken from Ref.P]. The energies found are 
-47.03MeV and -50.6MeV respectively. In Fig. 2 we 
compare the neutron density profiles of the ^n(0~^) drop 
obtained from AFDMC calculations with AU8' and AU6' 
interactions, and from self consistent mean field calcula- 
tions using the Skyrme 1' model, and the SKM model. 



O.lr 



0.08 -. 



Skyrme I' 
SKM 



-0.06--, 




FIG. 2: Density profiles for a *n drop for Vo = 20MeV. Solid 
line: result with the interaction; dashed line: result with 
the v^' interaction; dotted-dashed line: HF with Skyrme 1'; 
dotted line: SKM. 



As it can be seen, both of the energy density functional 
results differ significantly from the quantum Monte Carlo 
ones, particularly in the peak density value, at the center 
of the drop. The SKM however shows a density which is 
much closer to the quantum Monte Carlo result than the 
Skyrme 1' one. 

In Table ITTll we report the results of AFDMC calcula- 
tions for a *n droplet for different values of the depth of 
the confining potential Vq. The importance function for 
each case was optimized only with respect to the param- 
eters of the potential well from which the orbitals to be 
included in the Slater determinant are obtained. 

Together with the AFDMC energies we report the vari- 
ational energies E^- The difference between tha varia- 
tional and the AFDMC energy is quite large, in all cases 
between 20 and 30%. This indicates that the evaluation 
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TABLE III: AFDMC ground state energy (in MeV) for a '^n 
droplet as a function of the parameter Vo (in MeV) of Vext- 
The energy difference AE = E{AU8') ~ E{AU6') gives an 
unbiased estimate of the spin-orbit interaction contribution. 



Vo 


Et 


Eaus' 


{V) 


{'"l-s) 


Eau& 


A£ 


20 


-26.27(7) 


-37.8(2) 


-124.6(3) 


-0.41(7) 


-36.82(5) 


-1.0(2) 


25 


-48.0(3) 


-60.5(1) 


-168.2(4) 


-0.59(4) 


-59.6(1) 


-0.9(2) 


30 


-69.3(3) 


-85.3(2) 


-201.3(2) 


-0.72(4) 


-83.7(1) 


-1.6(4) 


35 


-92.7(3) 


-111.0(2) 


-232.4(2) 


-0.71(3) 


-109.86(5) 


-1.1(4) 



0.015 



0.010 



Q- 0.005 - 



0.000, 




of mixed estimators is at the same level of accuracy for 
all the cases, and the observed trends should not be too 
affected from the particular choice of the trial functions. 
The mixed energies are given for the AU9>' and AUQ' in- 
teractions, in order to calculate the contributions of the 
spin orbit interaction to the total energies. These are 
given in by /S.E on the last column of the figure. 

We have also evaluated the mixed estimate of the spin- 
orbit interaction. Since {^\vi g\^) = 0, the spin-orbit 
expectation value (v^.g) is approximately given by half 
the mixed estimate, and is reported in the fifth column 
of the Table. It is not very sensible to the depth of the 
potential well, in contrast with the total energy, which 
varies almost linearly with it. 

One can see that (w^.^) is about 60% of AE in the 
range of depths of the well considered. On the other hand 
this value is very close to the spin-orbit splitting evalu- 
ated in the ''n drop reported in Table II. It should also be 
noticed that in the latter case, the values of the spin-orbit 
splitting obtained with the AU8' and the AU6' are very 
similar. This means that the evaluation of the spin-orbit 
splitting in is not strongly influenced by the spin- 
orbit term in the interaction. The spin-orbit splitting 
found in our AFDMC is smaller (about one half) than 
the values reported from GFMC calculations. However, 
the consistencies in our results mentioned above makes 
us confident in the robustness of the evaluation of spin- 
orbit contributions in our method. 

In order to understand the connection between the 
parameter Vq and the neutron density inside the drop, 
we have plotted in Fig. 3 the mixed estimators for 
the density profiles for all the confinements considered. 
Although the density at the center remains in a range 
0.06 < p{Q) < 0.09fm~'^, it shows a larger dip at the cen- 
ter, and, consequently a stronger localization of neutrons 
at a distance of ~ Ifm, for larger values of Vq. 

V. CONCLUSIONS AND PERSPECTIVES 

In this paper we have applied the Auxiliary Field Dif- 
fusion Monte Carlo method to the study of a closed-shell 
neutron drop ^n(0+) and of the open shells '^n{l/2^) and 



FIG. 3: Neutron density profiles in a n drop for different 
values of the depth of the confining potential Vo 

^n(3/2^) confined by an external Wood-Saxon field. The 
method compares rather well with standard GFMC tech- 
niques. Since the AFDMC, scales better than GFMC 
with the number of particles it is currently the only vi- 
able quantum Monte Carlo method to study large neu- 
tron rich systems. 

We have made a detailed analysis of the contribution to 
the total energy of the spin-orbit interaction, which pre- 
vious calculations on neutron matter [l^ ha-s indicated 
as the main difference between AFDMC and GFMC cal- 
culations. Such analysis has been carried on at various 
values of the depth Vq of the confining potential V^xty 
and, therefore, for different confined neutron structures. 
We have found a spin-orbit contribution which is gener- 
ally smaller than in GFMC calculations, confirming the 
findings of ref. [ifij . Such contribution lowers the binding 
energy of about 1 MeV in all considered cases. 

The results obtained for the binding energies of the 
■^71(1/2-) and "^71(3/2") droplets confirm earlier GFMC 
results on the spin orbit splittings, which are significantly 
smaller than those provided by the energy density func- 
tional models. The AFDMC value of the spin-orbit split- 
ting is about one half of the value predicted by GFMC. 

Future work includes improving the trial function VP, 
as well as the guiding function "^q for transient estimate. 
The properties of oxygen and flourine istopes can be cal- 
culated with the AFDMC method. 
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